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1. INTRODUCTION 


This report summarizes research supported at Oceanweather Inc. and carried out nominally over a 
two-year period beginning in early 1999 and comprising a collaborative effort involving research 
teams at three institutions: Dr. Vincent J. Cardone and co-investigator Andrew T. Cox from 
Oceanweather Inc., Dr. W. Linwood Jones and co-investigators at the University of Central Florida 
(UCF), and Dr. Willard J. Pierson and co-investigators at the City College of City University of New 
York (CUNY). The research addresses two general problems: (1) the demonstration and improvement 
of scatterometer surface wind retrievals in tropical cyclones; (2) diagnosis of tropical cyclone 
planetary boundary layer (PBL) wind fields using a combination of high resolution NASA 
Scatterometer (NSC AT and its successor QuickScat) retrievals and a primitive equation vortex PBL 
numerical model. This research addresses the objectives of the US Weather Research Program 
(USWRP) related to improving hurricane forecasts near landfall including nowcasting storm intensity 
and tendency and prediction of hazards such as extreme winds and storm surge in coastal areas. The 
ultimate goal is to provide operational facilities, such as the NOAA Tropical Prediction Center, with 
improved methods for analyzing hurricanes and forecasting their movement and impact. 

This report focuses on the contributions of Oceanweather Inc. to the collaborative program. Later 
this year an integrated team report will be submitted upon the termination of the part of program 
supported at CUNY. 

The specific problems that must be solved to achieve the long-range goal may be stated as follows. 

1 . Demonstrate that high-resolution scatterometer measurements over tropical storms and other high- 
marine surface wind regimes when processed through an appropriate geophysical model function 
possess sufficient accuracy and dynamic range to be useful indicators of cyclone intensity and surface 
wind structure. 

2. Develop a method which utilizes the high resolution scatterometer (from NSCAT and is successor 
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QuickScat) data acquired in a single pass over a tropical cyclone (TC) in a consistent way to diagnose 
basic storm intensity and size in terms of parameters routinely used by warning centers including 
maximum sustained surface wind, V mx , and the radii of 35 and 50 knots sustained winds, R 35 and 

^ 50 - 


3. Utilize scatterometer data acquired in a succession of passes over a TC in a consistent way to 
specify the complete time and space evolution of the surface wind field over the sea suitable to drive 
ocean response and wind damage loss models thereby improving predictions of hazards to life and 
property from extreme winds, storm surge and coastal surf. 

These objectives have been substantially demonstrated in the research completed to date as 
summarized in the next sections and documented in detail in publications and presentations supported 
by this program and additional papers in preparation. In addition, an important first step has been 
taken to implement the results of this research operationally by integrating the analysis methods 
developed in this study into a major team program initiated in 1999 under the Navy Ocean 
Partnership Program (NOPP) organized by Prof. H. C. Graber at the University of Miami. This 
NOPP program is entitled: “Real Time Forecasting of Winds, Waves, and Surges in Tropical 
Cyclone”. The program includes the endorsement and participation of the NOAA Tropical Prediction 
Center and the NOAA/AOML Hurricane Research Division. Oceanweather Inc. are a critical 
member of the NOPP and are primarily responsible for developing the module to be used for 
operational analysis and prediction of the tropical cyclone wind fields. The early results of this 
program (a four-year NOPP program is planned) are to be reported in a paper in preparation for the 
AMS Annual Meeting, January 2001 (Abstract attached as Appendix D). Finally, we should note that 
the success of the early results of our research has evidently stimulated significant interest in the 
scientific community on the utilization of scatterometer data to study tropical cyclones. 
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2. DEVELOPMENT AND TESTING OF A NEW SCATTEROMETER GEOPHYSICAL 
MODEL FUNCTION AND WIND RETRIEVAL METHOD 


2.1 Analysis of NSCAT in Typhoon Violet. 

The model function development work was carried out mainly by Professor Jones and his colleagues 
after it was found that the NSCAT Project Baseline geophysical model function, NSCAT1, yielded 
wind speeds in tropical cyclones which were too low. As part of an early test of the new function 
Oceanweather developed a wind field for typhoon Violet in three NSCAT passes to develop the 
comparison database for the NSCAT evaluations. The wind fields utilized all conventional data 
available for those cases and its PBL vortex model. The new model function was developed from 
aircraft scatterometer data in a tropical cyclone and other supplementary data that indicated that the 
backscatter values measured in a tropical cyclone are much lower for a given wind speed and 
direction than the backscatter values were used in NSCAT- 1 . The new model function, called a 
neural net geophysical model function, was found to yield improved wind speeds in Typhoon Violet 
when applied with a Spatial Adaptive Retrieval Algorithm (SARA) which uses a-priori wind 
directions derived from Oceanweather’ s PBL model analysis, and a goodness of fit algorithm to filter 
out possible rain contaminated cells. SARA allows each individual backscatter value to be used to 
retrieve a wind speed. 

The analysis of the validation wind field in Typhoon Violet and the SARA algorithm evaluation is 
documented in the following paper which appeared in the Journal of Geophysical Research special 
issue: 

Jones, W.L., V. J. Cardone, W.J. Pierson, J. Zee, L.P.Rice, A.T. Cox, W.B. Sylvester, 1999: 
NSCAT High Resolution Surface Winds Measurements in Typhoon Violet. J. Geophys. Res ., 
NSCAT Special Issue, 104, C5, 11247-11259. 

This paper presented the analysis of NSCAT data in Revs. 478 and 485 which sampled much of the 
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circulation of Western North Pacific Typhoon Violet, a tropical cyclone that earlier had attained 
super-typhoon intensity but that at the time of these revs was of lower intensity with maximum 
sustained wind speeds of 90 knots. In the absence of aircraft reconnaissance data, the analysis of the 
surface wind “ground truth” relied on conventional data and the PBL model. The only use of 
NSCAT data in the surface wind analysis was a slight repositioning of the eye based on analysis of 
NSCAT sigma-0 profiles along constant incident angles by the CUNY group and research carried out 
by Mladen Susanj of the UCF team as described in more detail in the UCF report on this project. 

The NSCAT data compared in Violet benefited from three new analysis features. First, the 
geophysical model function (called NMGMF at the time) was refitted and improved over the standard 
NSCAT- 1 function by using backscatter measurements and coincident aircraft measurements of wind 
in hurricanes reduced to the surface using empirical flight-level 10 m reduction factors (Donnelly et 
al., 1998). According to the new function, the sigma-0 backscatter values are lower than NSCAT- 1 
for wind speeds greater than about 25 m/s but they increase very slowly with increasing wind speed 
even to wind speeds up to 40 m/s. 

The second novel analysis feature is the Spatial Adaptive Retrieval Algorithm (SARA). Whereas the 
normal NSCAT processing produces wind vectors located on 50 km centers from groups of NSCAT 
cells, SARA recovers wind vectors at the intrinsic resolution of each sigma-0 cell, namely 8 km x 35 
km cells sampled on 25 km centers along the satellite track. The wind direction alias removal step is 
averted because the field of vector wind direction is tightly constrained about a tropical cyclone and 
may be specified from a first guess PBL model solution as solved from conventional data only. 

The third feature was the application of a “goodness of fit” algorithm to flag retrievals probably 
contaminated by heavy rain. Basically, this algorithm examines the anisotropy of four nearest 
neighbors of sigma-0 to determine whether one or more of the measurements are inconsistent within a 
25 km grouping because of rain. This algorithm is described by Jones et al. (1999). The application 
of this algorithm to Violet indicated that even close to the main inner core of the circulation there are 
cells successfully retrieved because the beam is able to penetrate to the surface in areas of relatively 
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light rain or clouds between the spiral rain bands. 


Figures 1 shows the retrieved wind field plots and Figure 2 shows scatter plot comparisons of 
retrieved and analyzed surface winds in the two Violet revs with the best coverage of NSCAT data. 
The field plots indicate that the SARA algorithm provides significant improvement in spatial 
sampling over the standard NSCAT processing over the high wind areas near the center. The scatter 
plots show that at least qualitatively, the PBL model winds are in good agreement with the new 
NSCAT retrievals with winds recovered with little bias up to wind speeds to 35 m/s although the 
NSCAT wind speeds still slightly underestimate the model winds above about 30 m/s. These plots 
refer to wind speeds as the “average” wind at 10 meters. Operational warning centers classify 
tropical cyclones in terms of “sustained” wind speeds, which are about 25% greater than the average 
wind speed. Our study of NSCAT in Violet demonstrates, therefore, that the “sustained” winds up to 
perhaps 45 m/s or about 90 knots may be sensed even close to the center. 

2.2 Analysis of Hurricane Lili 

The analysis of the PBL wind field in Typhoon Violet for the selected passes in October, 1996 could 
not benefit from the best in-situ wind data, namely reconnaissance aircraft data, because since 1986 
there has been no aircraft reconnaissance of North Pacific Ocean tropical cyclones. During the 
formulation of this project in the Spring of 1997, it was expected that NSCAT would acquire data in 
several North Atlantic hurricanes during the 1997 season coincident with aircraft data. However, the 
early demise of NSCAT on June 30, 1997 precluded this prospect. Therefore, a thorough search was 
made for potentially useful NSCAT data acquired over tropical cyclones during its lifetime. 

Oceanweather carried out this search using all available meteorological data and the outputs of 
various warning centers while UCF scanned the NSCAT data set. We were pleasantly surprised 
during this search to find a number of potentially interesting cases in the North Atlantic Ocean, 
including Gulf of Mexico and Caribbean Sea, including Tropical Cyclone Josephine, October 4-9, 
1996 and Hurricane Lili, October 15-21, 1996. Because Tropical Cyclone Josephine exhibited 
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properties of extratropical structure, Hurricane Lili was selected for further analysis. 

Hurricane Lili developed in the western Caribbean and moved across Cuba and into the southwest 
North Atlantic Ocean. Lili reached Category 3 intensity while under aircraft surveillance south of 
Bermuda with central pressure of 960 mb and maximum flight level wind speeds of 1 1 0 knots. 

A search of NSCAT data showed that Hurricane Lili was an excellent case for analysis and it was 
agreed by the team that Oceanweather should proceed to an intensive analysis of the PBL winds in 
this storm at the times of Revs 900 (0342 UTC, October 19), Rev 907 (1432 UTC October 19) and 
Rev 914 (0314 UTC October 20). The best-documented rev was rev 900 on October 19, 1996 
because of the wealth of aircraft data. Figure 3 (a-d) shows the coverage of aircraft data within +/-6 
hours of this rev and gives a comparison of the aircraft data and modeled winds for four different 
flight level- 10 m reduction algorithms. The analysis of NSCAT data in this rev is presented in the 

following paper: 

Cardone, V.J., A.T. Cox, W. J. Pierson, W. B. Sylvester, W.L. Jones and J. Zee.: NSCAT 
Scatterometer High Resolution Wind Fields for Hurricane Lili. IEEE IGARRS’99, Hamburg, 
Germany, June 28-July 2, 1999, 1509-1513. 

Figure 4 (from this paper) shows the modeled wind field solution and the retrieved winds using a 
new Tropical Cyclone Geophysical Model Function (TC-GMF) (Zee and Jones, 1999) and the SARA 
algorithm. Figure 5 shows the comparison of the retrieved winds speeds and the model wind speeds 
at each SARA cell which passed the rain filter. This comparison confirms the basic properties of the 
NSCAT winds with the new model function, namely good wind speed sensitivity up to wind speeds 
of at least 20 m/s and a small bias above 25 m/s. Above 25 m/s, the bias in the retrieved wind speeds 
is only 0.1 m/s or TC-GMF compared to a bias of -5. 32 m/s or NSCAT-lb. 
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2.3 North Sea NSCAT vs. Offshore Platform High Wind Speeds 


The new model function used for the analyses summarized above was derived from the assumption 
that the NSCAT- 1 model function is correct in the low to moderate wind speed range and used 
aircraft scatterometer data sets at higher wind speeds (up to 35 m/s). Beyond 35 m/s the new 1 
function is an extrapolation. An important difficulty is this approach is that the aircraft winds are 
reduced to 10 m from flight level using an empirical algorithm. This process introduces an 
uncertainty and potential bias in the reference 10 m wind speeds above about 20 m/s. Meanwhile, the 
validation of scatterometer winds have generally utilized comparison data sets assembled from data 
buoys and research vessels that are generally limited to 10- m neutral wind speeds of 20 m/s or less. 
For example, Freilich and Dunbar (1999) found only 184 out of 56,000 NSCAT-buoy collocations 
with 10-m buoy winds exceeding 1 8 m/s and virtually none above 22 m/s. Validation is further 
complicated by evidence that buoy wind speeds may be biased low in high sea states. Our evaluation 
of the sensitivity and accuracy of scatterometer winds above 20 m/s or so by comparison of winds 
from the cyclone PBL model in tropical cyclones is encouraging; however, the lack of high- wind 
speed in-situ data sets coincident with NSCAT and QuickScat hampers further refinements of the 
backscatter model function. 

To ameliorate this lack of in-situ data at very high wind speeds, as part of this research program we 
assembled a unique comparison database for NSCAT wind validation at high wind speeds. This data 
set utilized winds measured from offshore platforms in the open Northern North Sea and Southern 
Norwegian Sea. Of the dozens of platforms which make and report synoptic observations in this 
region, we selected only nine platforms at each of which a calibrated anemometer is well exposed at 
the top of the drilling derrick (the height range is 86-m to 143-m) and quality-controlled continuous 
time average wind samples are archived by and available from European weather centers. 

Collocation of the NSCAT 25-km database with the platform database yielded 3663 (2773) matches 
for separation distances of 25-km and time offsets up to 1.5 (0.5) hours. Using Liu's reduction to 10- 
m equivalent neutral wind, 1 10 matches are found for wind speeds above 20 m/s and 37 matches for 
wind speed above 24 m/s. However, because in some regimes the anemometer may be above the 
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surface boundary layer, four different stability-dependent wind profile models were used in the 
analysis, including a model that extends the wind profiles into the Ekman layer. The data were 
analyzed and reported to the AGU 2000 Spring Meeting, May 30-June 3, 2000. Washington, DC 
(Cardone et al., 2000). 

The NSCAT-platform winds were compared in terms of scatter plots, regressions on bin-averaged 
data and comparisons of wind speed and platform wind speed distributions in terms of quantile- 
quantile scatter plots. For a selected set of 58 passes during particularly intense storms, the number 
of collocations was greatly expanded (about 10,000 matches) by analyzing the platform data into a 
continuous field using the natural neighbor objective analysis system. Figure 6 shows one of these 
high- wind speed revs (February 16, 1997, 2100 GMT) and Figure 7 shows the time history of 
measured winds at one site (drill-ship “Polar Pioneer”). Average wind speeds up to 45 m/s were 
measured at anemometer height of 92 m in this event and were reduced to 1 0 m using the four 
different profile forms noted above. Figure 8 shows the NSCAT (using NSCAT-1 model function) 
versus platform wind speeds over all available comparisons (all revs in the mission). The standard 
deviation of the difference in wind speeds is only 1 .5 m/s but there is a clear low bias in the NSCAT- 
1 winds above about 15 m/s. This high wind speed negative bias is also shown in the distributional 
comparisons (quantile-quantile plots) given in Figure 9 for the difference profile forms. This pilot 
study confirms that North Sea platform data may be used to supplement cyclone aircraft data to refine 
the model function at high wind speeds. Unfortunately, the four different profile forms used do not 
yield the same 10-m wind speeds for a given stability and anemometer level measurement. 
Experiments are underway in Canada to resolve the difference between platform and buoy derived 
estimates of the 10-m wind and in the UK to study the flow distortion (if any) contribution to winds 
measured from offshore platforms. 

During the Northern Hemisphere winter sensed by QuickScat, 1999/2000, a number of very intense 
and damaging storms crossed from the North Atlantic Ocean over the North Sea and into NW Europe 
and it is therefore possible to produce a data set such as assembled here for validation of QuickScat 
and further refinement of the model function. 
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Appendix A gives copies of key viewgraphs presented the AGU conference. A paper is being 
prepared for submission to the next JGR special issue as follows: 

Accuracy of Scatterometer Winds Assessed from In-Situ Measured Wind Data Up to 32 m/s 
VJ.Cardone, E.A. Ceccacci, A.T. Cox, J. G. Greenwood 



3. INVERSE MODELING OF TROPICAL CYCLONE WIND FIELD FROM QUICKSCAT 


3.1 Approach 

The inverse model developed and tested in this study to diagnose the inner core structure of a tropical 
cyclone from scatterometer data utilizes the most recent update of Oceanweather’s tropical cyclone 
PBL model (Thompson and Cardone, 1 996) hereinafter TC96. This model solves the primitive 
equations of motion representing the vertically integrated PBL on a nested rectangular grid in a 
moving coordinate system centered on the storm. Vertical and lateral friction are included. The model 
resolution of the inner core structure on the inner grid nest is 2 km . The model is driven by a 
specification of the following parameters (see TC96, p. 202 for a theoretical description of the model 
and its initialization): 

speed and direction of vortex motion 

equivalent geostrophic flow of the ambient PBL pressure field 
in which the vortex propagates 
total storm pressure anomaly 

pressure anomaly associated with the first component of 
the exponential radial pressure profile 

scale radii of the up to two components of the exponential radial pressure 

Holland’s (1980) profile peakedness parameter for each component 

The problem may be posed as follows. Given only scatterometer winds acquired in a pass over the 
TC, standard warning center satellite derived Dvorak estimates of intensity (as a first guess) at the 
time of the pass, and conventional synoptic data available on the periphery of the circulation, how 
may TC96 be used in an "inverse" sense to develop refined estimates of storm intensity (i.e. 
maximum sustained eye-wall wind speeds) and size, namely R^, Rp 2 R35 , R 50 By inverse modeling, 
we mean an objective algorithm that seeks and selects the optimum TC96 forcing variables above. 


Yjpj, Vj ]r 

V gs , D g s 

Dp 

dpi 

RpU Rp2 

profile 

B], B 2 
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which when applied to TC96, return a surface wind field that exhibits minimum difference, in a least 
squares sense, between the model winds and the scatterometer winds. Such an inverse modeling 
procedure has already been constructed (Car done et al. 1994) for the case when only dropsonde 
central pressure and the radial profile of the flight level wind as measured by an aircraft are known. 
That problem is somewhat simpler than the scatterometer problem because in the scatterometer case 
the maximum eye-wall wind speed is not known if the wind speed is in the model function saturation 
range. We assume here that rain contamination and possibly model function saturation for the more 
intense storms prevents direct retrieval of this quantity (it is approximated by the Dvorak estimate), 
whereas when there is aircraft reconnaissance of a cyclone, the aircraft almost always penetrates the 
eye-wall on a radial leg and accurately measures the flight level wind in the eye-wall. 

The preliminary results carried out with the neural net algorithm summarized above suggests that 
average wind speeds up to 25 m/s may be routinely retrieved with little bias in tropical cyclones, but 
it is not yet clear that wind speeds above 30 m/s will be retrievable. Therefore, we have developed a 
preliminary design for an inverse model which seeks to find solutions of TC96 which best match the 
azimuthally averaged scatterometer wind field in areas of wind speeds of less than 20 m/s (1-hour 
average at 1 0 m, not sustained). This part of the circulation includes the areas of the operationally 
significant radius of the 35 knot and 50 knot sustained winds. The inverse model required the 
generation of over 30,000 TC96 model solutions covering the following model parameter (single 
exponential only) space: 


Latitude 

^spd’ ^dir 

V S s, D gs 

D P 

Rr 

B 

H 

Dif 


20 degrees 

stationary 

stationary 

2 to 120 mb in 2 mb increments 
5 nm to 1 20 nm in 4 nm increments 
1 to 2.5 in 0.1 increments 
500 m (boundary layer depth) 

-2 C (air-sea temperature) 
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For each solution the following information was saved: 


Maximum surface wind 
Radius of maximum wind 
Azimuthally averaged wind speed in 5 nm bins 
Radius (maximum) of 35 knot winds 
Radius (maximum) of 50 knot winds 

A preliminary working algorithm was applied which matched only scatterometer and model radii of 
35 and 50 knot winds but we found a better solution was obtained by matching the entire azimuthally 
averaged wind profile over the valid range of wind speeds in an rms sense. 

3.2 Application of Inverse Model to Floyd (1999) using QuickScat 

3.2.1 General History of Floyd 

Floyd was a large and very intense Cape Verde type hurricane that formed on September 7, 1999 in 
the eastern North Atlantic and later seriously threatened the southeast US coast triggering evacuations 
of millions of people. Floyd attained the high-end of Category 4 status with peak sustained winds of 
135 knots east of the Bahamas on September 13, 1999. The system weakened slowly thereafter and 
induced immense flooding as it entered North Carolina on September 1 6 with peak sustained wind 
speeds of about 90 knots. Floyd caused 56 deaths and 3-6 billion dollars of property damage in the 
United States. 

Almost continuous aircraft monitoring of Floyd from September 8 through the 16 th defined the 
position, central pressure, maximum eye-wall wind speed and two-dimensional wind field structure 
of the storm. Those data allowed the NOAA AOML Hurricane Research Division to construct 31 
“snapshots” between September 9 and 16 of the sustained wind field at 10 m over water using a 
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method (Powell et al., 1998) which consists of objective analysis of mainly surface winds 
extrapolated from the aircraft flight level winds using empirical reduction factors. Figure 10 for 1330 
UTC 13 September gives an example of these wind fields. We have used these snapshots to construct 
a time and space continuous representation of the wind field to serve as a basis for the evaluation of 
wind fields inverse modeled from the QuickScat data. 

3.2.2 Inverse Model Fits 

All QuickScat passes over the North Atlantic between September 8 and 16 were downloaded from the 
JPL distribution center. It was found that because of the wide swath of QuickScat, at least part of the 
storm circulation was viewed at least twice per day, with pass times typically near 00, 09, 12 and 21 
UTC. Appendix B shows the coverage of the QuickScat data for each of the 32 data swaths for which 
there was data within about 360 nm of the center of the storm. For about half of these passes, the 
center of Floyd and its entire inner core circulation are well sampled except, of course, for the gaps 
which indicate the rain flag was set (we used the preliminary QuickScat processing for this study, not 
the more recently reprocessed data which only became available during summer of 2000). 

The QuickScat data were azimuthally averaged within 360 nm of the center in 5 nm bins (we used 
the official fix, though in a case without aircraft data, satellite images or the sophisticated processing 
of the scatterometer data developed at UCF may be used to infer a center). Since we are using 
standard processed QuickScat data that uses essentially the old NSCAT-1 model function there is the 
possibility of attenuation/saturation of the winds speeds above 20 m/s. Only QuickScat bins with 
azimuthally averaged wind speed less than 20 m/s were used as a conservative way to omit 
QuickScat data which might have saturated. Also, passes for which no data were retrieved within 
150 nm of the center were discarded. Only 7 of the passes were discarded for this reason. 

As noted above, the simple inverse model developed involves matching the azimuthally averaged 
QuickScat radial wind speed profile with precomputed TC96 model azimuthally averaged radial wind 
speed profiles in a least squares minimization sense given only the location of the center and the 
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equivalent Dvorak estimate of eye-pressure. Figure 1 1 shows the best match near 0000 September 
10, early in the life of Floyd when wind speeds even near the center were less than 20 m/s (hourly 
average). The profile radius and shape parameters selected are indicated on this plot. Figure 12 is an 
example match near 2100 UTC September 1 1 . In this case even though the azimuthally averaged 
QuickScat winds greater than 20 m/s were not used, the best- matched profile gives azimuthally 
average peak winds of about 29 m/s at a radius of maximum wind (Rmx) of about 35 nm. 
Interestingly, the QuickScat data in the vicinity of Rmx also average to about 30 m/s. However, the 
case for Floyd near peak intensity shown in Figure 13 indicates that the inverse model provides 
maximum azimuthally average wind speeds of 47 m/s while the QuickScat wind speeds there average 
only to about 30 m/s. In fact, it would appear that within 120 nm the of the center, the QuickScat 
winds in this pass are biased low but there is sufficient information in the profile outside that radius to 
allow the inverse model to recover the inner-core structure. All of the inverse model fits developed 
are included in Appendix B. 

Each of the inverse model fits provides the critical model “snapshot” profile parameters. To these 
must added the storm motion vector (taken from the past track) and the background pressure field 
(taken from standard large scale NOAA NWP fields). Then the TC96 model is solved for each 
snapshot of the actual wind fields and these snapshots are interpolated to a regular grid of spacing of 
0.125 degree and to 30-minute time intervals using OWI’s standard software for this purpose. 

3.2.3 Comparison of Inverse Model Wind Fields and QuickScat Winds 

The snapshots wind fields corresponding to each QuickScat pass may be compared to the individual 
QuickScat retrieved winds. This direct comparison is shown in Figure 14 for all snapshots recovered. 
As expected, the higher QuickScat winds are lower than the model winds. Over these 13,502 
comparisons, the mean difference (Quickscat inverse model field - QuickScat individual retrievals) is 
0.81 m/s and the standard deviation of the difference is 3.63 m/s and much of this scatter arises in 
data pairs above 20 m/s in which range the QuickScat data are biased and relatively inaccurate. 
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3.2.4 Comparison of Inverse Model Wind Fields and HRD Wind Fields 


While the HRD winds do not constitute an absolute standard for verification, they do provide a 
measure of the success of our inverse model to utilize QuickScat data to diagnose the complete two- 
dimensional surface wind field in a severe tropical cyclone. Comparisons of HRD derived wind 
fields and inverse model wind fields are presented at the times of the 24 HRD snapshots in Appendix 
C between September 1 1 and September 1 6. Each wind is shown as a vector and wind speed 
differences are color-coded. The statistical differences between the HRD and inverse model winds 
speeds and directions over the domains shown and for the times shown in Appendix C are given in 
Table 1 . There are generally between 6000 and 7000 grid point comparisons for each snapshot. For 
example, near the peak of the storm on September 1 3 there are four snapshot comparisons at around 
00UTC, 13UTC, 19UTC and 22 UTC. The mean differences in wind speed (Inverse model-HRD) 
are 3.23 m/s, -.2.3 1 m/s, 0.61 m/s and 1 .00 m/s respectively. The standard deviations of the wind 
speed differences are: 2.34 m/s, 4.38 m/s, 2.70 m/s and 2.57 m/s while the correlation coefficients for 
wind speed are .94, .80, .91, and .92. The mean wind direction differences are -.75 deg, -3.29 deg, - 
8.36, -1 1 .42. These differences are in the sense that there is more inflow in the model winds than in 
the reduced aircraft winds and this may well reflect insufficient turning of the aircraft winds toward 
the surface in the HRD analyses. Over all snapshots, the differences are: 

Wind speed bias (QuickScat inverse model-HRD) 2.79 m/s 


Wind speed difference standard deviation 2.63 m/s 

Wind speed scatter index . 1 7 

Wind speed correlation coefficient .92 

Wind direction mean difference -5.6 degrees 

Wind direction difference standard deviation 1 6. Degrees 


These statistics indicate a close statistical match between the QuickScat derived inverse model winds 
and the HRD aircraft derived winds. The difference plots of Appendix C indicate that the wind 
differences within a given snapshot display occasionally large spatially coherent differences, and 
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these may arise from a number of causes including: (1) failure of the TC96 model to simulate smaller 
scale spiral bands which are indeed sensed by the aircraft; (2) HRD wind errors in storm quadrants 
which might not have been probed by the aircraft in the composite of data used for a given snapshot; 
(3) slight positioning differences of wind field features which have been qualitatively well modeled 
in both field. What is encouraging, however, is that the coherent spatial differences appear not to 
have much coherency in time from snapshot to snapshot, suggesting that such errors will have 
minimal impact on ocean response models forced by the QuickScat winds. 
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4. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER RESEARCH 


We may now restate the specific problems addressed and summarize the results achieved. 

Demonstrate that high-resolution scatterometer measurements over tropical storms and other 
high-marine surface wind regimes when processed through an appropriate geophysical model 
function possess sufficient accuracy and dynamic range to be useful indicators of cyclone 
intensity and surface wind structure. 

The analysis of NSCAT data sets in Typhoon Violet and Hurricane Lily have demonstrated that 
accurate scatterometer measurements of surface wind are possible even near the inner core of intense 
tropical cyclones for wind speeds up to about 35 m/s, but this dynamic range requires a geophysical 
model function different in its tuning at high winds than used operationally for NSCAT and 
QuickScat, and a spatial adaptive retrieval algorithm which utilizes the fact that the vector wind 
direction is so tightly constrained in the inner core of mature tropical cyclones that wind direction can 
be assumed to be prescribed accurately independently of the remote sensing data and utilized to 
resolve surface winds at the highest spatial resolution allowed by the scatterometer cell footprint. The 
filtering of rain-contaminated cells remains a daunting problem but the UCF goodness of fit 
algorithm used in this study allows useful high wind retrievals even near the cyclone center. 

However, operational NSCAT and QuickScat winds remain biased low at wind speeds above about 
20 m/s and we recommend additional studies to develop a more accurate model function at high wind 
speeds. A pilot study using North Sea platforms to validate NSCAT data during the winter 1996/97 
suggests that such data may be used to refine the QuickScat model function. 

Develop a method which utilizes the high resolution scatterometer (from NSCAT and is successor 
QuickScat) data acquired in a single pass over a tropical cyclone (TC) in a consistent way to 
diagnose the surface wind field in the inner core and far-field. 
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The coupling of a highly tuned and proven dynamical model of the cyclone PBL (TC96) with the 
operational QuickScat data has led to an inverse model which recovers a snapshot of the asymmetric 
surface wind field from a QuickScat pass over a cyclone provided the circulation within 150 km of 
the center is included in the swath. The wind field including the eyewall region is recovered with 
wind speed bias and scatter of about 2.5 m/s and wind direction bias and scatter of about 5 Degrees 
and 15 Degrees respectively, all relative to a reference wind field snapshot derived objectively from 
copious aircraft flight level winds reduced to 10 m.. The wind field allows accurate recovery of the 
storm intensity and size parameters routinely used by warning centers including maximum sustained 
surface wind, V mx , and radii of 35 and 50 knots sustained winds, R 35 and R 50 . In addition to the 
QuickScat data, the method requires knowledge of the position of eye of the storm, and an accurate 
Dvorak type estimate of intensity from which may be specified an estimate of eye pressure. We 
suggest that reprocessing of QuickScat with the improved model function developed at UCF will 
yield an even greater dynamic wind speed range and even more accurate inverse model parameters. 
We recommend the analysis of additional tropical cyclones after the improved model function is 
implemented. 

Utilize scatterometer data acquired in a succession of passes over a TC in a consistent way to 
specify the complete time and space evolution of the surface wind field over the sea suitable to 
drive ocean response and wind damage loss models thereby improving predictions of hazards to 
life and property from extreme winds, storm surge and coastal surf. 

It has been demonstrated through our Hurricane Floyd case study that QuickScat data may be used to 
inverse model the complete life-cycle of the surface wind field in a typical very intense North 
Atlantic hurricane. We recommend study of additional cases of this type, including application of the 
derived wind fields to force ocean response models for surface waves, continental shelf currents and 
coastal storm surge and validation of the ocean response predictions against measurements. 

This study represents an important first step toward the operational implementation of the analysis 
methods developed in this study. The operational aspects of this work are being continued under a 
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major team program initiated in 1999 under the Navy Ocean Partnership Program (NOPP) entitled: 
“Real Time Forecasting of Winds, Waves, and Surges in Tropical Cyclones”. The program includes 
the endorsement and participation of the NOAA Tropical Prediction Center and the NOAA/AOML 
Hurricane Research Division. The early results of this program are to be reported in a paper in 
preparation for the AMS Annual Meeting, January 2001 . 
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Table 1. Differences in Hurricane Floyd (1999) between OW1 10 m average wind fields inverse 
modeled from QuickScat and NOAA HRD 10 m average wind fields derived from aircraft data 
QuickScat derived wind fields are interpolated in space and time to the grid and times o t e 
24 available HRD snapshots. 
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SARA max wind retrievals using NNGMF model function and TC96 wind directions 
for Violet revolution 478, showing (left) SARA max wind vectors and (right) wind speed differences 
between TC96 and SARA max . Open triangle denote < -10 m s -1 , open squares represent -10 
to -5 m s -1 , and open diamonds denote -5 to 0 m s“ J . Solid symbols correspond to positive 
differences. 



SARA max wind retrievals using NNGMF model function and TC96 wind directions 
for Violet revolution 485, showing (left) SARA max wind vectors and (right) wind speed differences 
etween TC96 and SARA max . Open triangles denote < -10 m s” 1 , open squares represent -10 to 
-5 ms , open diamonds denote -5 to 0 m s" 1 . Solid symbols correspond to positive differences. 
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Comparison of the wind speeds recovered 
by S ARA max with corresponding tropical cyclone model 
speeds (TC96) for revolution 478. Symbols are mean 
values for 5 m s -1 bins: asterisks denote NNGMF model 
function and minuses represent NSCAT 1 model func- 
tion. Vertical lines show the standard deviations of the 
SARA retrievals about the means. 



Comparison of the wind speeds recovered 
by SARA max with corresponding tropici cyclone model 
speeds (TC96) for revolution 485. Symbols are mean 
values for 5 ms' 1 bins: asterisks denote NNGMF model 
function and minuses represent NSCAT 1 model func- 
tion. Vertical lines show the standard deviations of the 
SARA retrievals about the means. 
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Measured Wind Speed (nrVs) 


NSCAT/LILI 1996 

OWI Model Winds vs. Flight Level Winds (adjusted) 

.76*FLW to SFC Reduction Used 

FLW +/-6 Hours from NSCAT Rev900 (October 19, 1996 4:00 GMT) 
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Figure 3a 



Measured Wind Speed (m/s) 


NSCAT/LILI 1996 

OWI Model Winds vs. Flight Level Winds (adjusted) 

GPS FLW to SFC Reduction Used 

FLW +/-6 Hours from NSCAT Rev900 (October 19, 1996 4:00 GMT) 



- 79 * - 78 * - 77 * - 76 * - 75 * - 79 * - 78 * - 77 * - 76 * - 75 * 



Grid Number Mean Mean 01 ff RMS Stnd Scat Corr 

Station Point of Pts Meas Hind (H-M) Error Dev Index Ratio Coeff 


438 19.62 20.53 0.91 3.62 3.50 0.18 0.63 0.81 

565 96.94 119.70 -16.92 N/A 54.99 0.15 N/A N/A 
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Figure 3b 



Measured Wind Speed (rrVs) 


NSCAT/LILI 1996 

OWI Model Winds vs. Flight Level Winds (adjusted) 

GPS/BOUNDY FLW to SFC Reduction Used 

FLW +/-6 Hours from NSCAT Rev900 (October 19, 1996 4:00 GMT) 
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Figure 3c 


Measured Wind Speed (m/s) 


NSCAT/LILI 1996 

OWI Model Winds vs. Flight Level Winds (adjusted) 

BOUNDY FLW to SFC Reduction Used 

FLW +/-6 Hours from NSCAT Rev900 (October 19, 1996 4:00 GMT) 
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Fig. 1 . Tropical cyclone PBL model surface winds from TC96 at 10 m elevation for Rev. 900. 
Wind barbs are shown at the corresponding NSCAT o" locations. 



Fig. 2. NSCAT retrieved winds from the Tropical Cyclone Geophysical Model Function for Rev. 900. 


Figure 4 





Fig. 3. Wind speed comparisons between the TC96 model and 
NSCAT retrieved winds for Rev. 900. The filled circles are for 
TCGMF and the tilted squares are for NSCAT- 1 b. 
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Unstable Stratification Regime 
Selected Time Series: 9702 Polar Pioneer 

Official" Power Law * Cardone (1969) (SBL) * Liu et al. (1979) (SBL) & Brown and Liu (1982) (PBL) 



(s/lu) sm 
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NSCAT WS (kt) 


PLATFORMS vs. NSCAT COLLOCATED - 9609-9706 

Using WindFN 



PLATFORMS WS (kt) 
Best Fit: Y = 0.820X + 1.680 
Total Points: 3669 
Mean X: 18.864 
Mean Y: 17.143 
Mean Diff: -1.721 
Root Mean Square: 3.537 
Standard Dev. : 3.091 
Scatter Index: 0.164 
Ratio: 0.246 
Correlation Coeff: 0.949 
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Quantile-Quantile Scatter Plots of Collocated 
NSCAT-North Sea Platform 10 m Wind Speed 
+ 1.105*x A 1 ' 77] 'I — -0.915 + 1. 21 7*x x 1 
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ATTENTION: HURRICANE SPECIALISTS 
Hurricane Floyd 1330 UTC 13 Sept. 1999 

Max. 1-min sustained surface winds (kt) for marine exposure 

Analysis based on US AFRES recon. 700 mb winds adj. to sfc: 0700- 1 154 z, 

3 GPS dropwindsondes MBL adj to st'c: 0750, 0927, and 1 1 15z 
1330 z position extrapolated from 1 121 z wind center fix assuming 280° @ 12 kt, mslp = 921 mb 



Figure 10 


file://H:\NOPPS\Floyd99\HRD\floyd_l 99909 1 3 1 330.gif 


12/28/99 




Range(Nmi) 




\D O 



Range(Nmi) 






1 



20 30 40 50 

Typhoon Model Wind Speed (m/s) 


Figure 14 


]?/ 31 / 9 < 






Appendix A. Abstract and key slides of AGU Spring Meeting, 2000 Paper: 

Accuracy of Scatterometer Winds Assessed from In-Situ Measured Wind Data up to 32 m/s 



Accuracy of Scatterometer Winds Assessed from In-Situ Measured Wind Data Up to 32 m/s 


V.J.Cardone, E.A. Ceccacci, A.T. Cox, J. G. Greenwood (Oceanweather Inc., Cos Cob, CT 06807, (203- 
661-3091) oceanwx@oceanweather.com, 


Evaluations of the accuracy of scatterometer winds have generally utilized comparison data sets assembled 
from data buoys and research vessels that are generally limited to 10- m neutral wind speeds of 20 m/s or 
less. For example, Freilich and Dunbar (1999) found only 184 out of 56,000 NSCAT-buoy collocations with 
10-m buoy winds exceeding 18 m/s and virtually none above 22 m/s. Validation is further complicated by 
evidence that buoy wind speeds may be biased low in high sea states. The scant evidence of the sensitivity 
and accuracy of scatterometer winds above 20 m/s or so has been provided to date only by comparison of 
winds from NWP systems or model generated wind fields in tropical cyclones. This lack of high-wind speed 
validation hampers scientific applications of existing data sets from NSCAT and QuickScat and further 
refinements of the backscatter model function. 


We have assembled a unique comparison data base for NSCAT wind validation at high wind speeds, utilizing 
winds measured from offshore platforms in the open Northern North Sea and Southern Norwegian Seas. Of 
the dozens of platforms which make and report synoptic observations in this region, we selected only nine 
platforms at which a calibrated anemometer is well exposed at the top of the drilling derrick (the height range 
is 86-m to 143-m) and quality-controlled continuous time average wind samples are archived by and available 
from European centers. Collocation of the NSCAT 25-km database with the platform database yielded 3663 
(2773) matches for separation distances of 25-km and time offsets up to 1.5 (0.5) hours. Using Liu's 
reduction to 10-m equivalent neutral wind, 110 matches are found for wind speeds above 20 m/s and 37 
matches for wind speed above 24 m/s. However, because in some regimes the anemometer may be above 
the surface boundary layer, four different stability-dependent wind profile models were used in the analysis. 
The data are analyzed and will be reported in terms of scatter plots, regressions on bin-averaged data and 
comparisons of wind speed and platform wind speed distributions in terms of quantile-quantile scatter plots. 


For a selected set of 58 passes during particularly intense storms, the number of collocations is greatly 
expanded (about 1 0,000 matches) by analyzing the platform data into a continuous field using the natural 
neighbor objective analysis system. These comparisons are reported as well. A similar analysis is planned 
using the reprocessed QuickScat following assembly of the platform data for the 99/00 winter. 
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though field experiments and research underway to resolve this issue 
(for a recent review see Taylor et al., 1999, CLIMAR99) 



PLATFORM WIND DATA 
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Sampling rate varies from 20 minutes (3 sites), 1-hour (3 sites), an 
hours (3 sites); averaging interval generally 10 min 

Owners: Shell, BP-Amoco, Phillips, Elf, Statoil, Norsk Hydro 









Platform Data Processing 
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Stable Stratification Regime 
Selected Time Series: 9706 Gullfaks 
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Collocation Process 
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Always match single nearest 25 km SWV within time and space filter 

Found nearly 3,700 matches: 30% in stable regimes, 70% neutral or 
unstable regimes 



NSCAT WS (kt) 


FRIGG vs. NSCAT COLLOCATED - 9609-9706 

Using WindBrown 
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PLATFORM WS (kt) 
Best Fit: Y = 0.909X + 0.245 
Total Points: 549 

Mean X: 18.342 
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Root Mean Square: 2.673 

Standard Dev. : 2.263 

Scatter Index: 0.123 
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PLATFORM WS (kt) 
Best Fit: Y = 0.769X + 2.308 
Total Points: 381 

Mean X : 19 • 840 

Mean Y: 17.556 

Mean Diff: -2.283 

Root Mean Square: 4.542 
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Scatter Index: 0.198 
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PLATFORMS vs. NSC AT COLLOCATED - 9609-9706 

Using WindFN 



Total Points: 3669 
Mean X: 18.864 
Mean Y : 17 . 143 
Mean Diff: -1.721 
Root Mean Square: 3.537 
Standard Dev.: 3.091 
Scatter Index: 0.164 
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Correlation Coeff: 0.949 
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Collocated NS CAT/North Sea Platforms 10 m Wind Speed Plots 

Binned on 1 m/s Pl atform Winds 
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Quantile-Quantile Scatter Plots of Collocated 








Mean Difference (Q-Q Regressions) 
NSC AT - Platform (10 m) 



Based on regression of 30 NDBC Buoys 
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All 58 Snapshots, WindFN adjusted with boundary 


Best Fit: Y = 0.900X + 

Total Points: 9992 
Mean X: 32.056 

Mean Y: 29.171 

Mean Diff (Y-X): -2.885 

Root Mean Square: 4.726 

Standard Deviation: 3.744 

Scatter Index: 0.117 

Ratio (Y>X): 0.191 

Correlation Coeff: 0.894 
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Natural Neighbor (nnGridR) WS vs. NSCAT 
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2. Resolution of wind profile 5 m - 150 m :platform-buoy 

comparison programs planned 

NSCAT sensitive to surface wind to at least 32 m/s with 
low bias of 5-15 % above 20 m/s 




Appendix B. QuickScat data coverage of Hurricane Floyd (1999) and best fits of inverse model 
parameters to azimuthally averaged wind speed (average wind speed at 10 m) 
































































Appendix C. Comparison of OWI winds from inverse model and HRD snapshot wind fields at 
time of HRD snapshots. Color contours of wind speed differences 




if) 00 







i <> ft: 


25 'QoCQB^i^ 0 !! QyVl Bi^LV^^r^s ^^ing^urricane 


- i 


/ y 
' / 
4 4 
4 4 
4 4 
4 4 
( ( 
I i 
l ( 
f i 
I i 


l l 



v 

V V V ■ B 


« 


■ 

■' ■' ■ X it. < B / 1 

V. 

v. V * B m 

a. 

rt 

« 


■ '. B". 'N/.*, SK at- It 1 Jf* J6 

v * 

\. v a ■. B 

a 

a 

a'. 


■' ■', --A *-t ■'• ■ ' *’- ■' ■ 

s. 

v A A, B 1. 




™ i 

■ , ■ ' if 

V 

V. V. V. B. ■. 




“ t 

■ , *, a 3 - a" ;/ 

V 

V V ■.. 





■-< ;*-*■»— B"* R' a* it" vC 

V 

V <5. V a.. 



b< 




V 

X. 


x. V X* 

V V V 
\> v 


;: i 

7 f 

f t 

‘ / 

/ / 

7 / 


■ ■■■■■■ ■ 
»»■■■ ■ ■ ■ 


v ^ ^ v v y y 

^ ^ ^ ✓ f s ■ / 

*t_ -±_ .>*_ /!_>/!_ /L /_ 


sne_F!oy^ 

1 QQQ 


" 1 1 
— A .A Jji.' 

■Q jj D> 

xtn 

.> 

> > 


Ji. 


A 

.A 

.A .A -A 

• A 




A 

*iw 

.A 

.A *\ 


- 

> 

> > 

.s* 

-- 


.A 

.A .A .A 

.A -A -% 

-* 

.■ 5 - 


> 

x 



.A . A .A 

-A -A -A 

• A 




x 


A 

»A .A .A 

-A -A * A 


-- 



-- 

X 

- 

.A .A -A 

-A *A -A 

.A 



- 

-- 


- 

.A -A .A 

• A *A A 

- 




- 

•- 

- 

- - - 

-X -X. *4 

-A 



- 

.w 

- 

- 

A A 

•A -A ‘V 




— 

-- 

*- 

*- 

- - - 

-x -x -x 

- ^ 

-- 

-- 

- 

- 

- 

- 

- - - 

*V *X ’X 

■A 




- 

- 

- 

-X 'X X 

v -v -v 

-A 



— 

— 

— 

- 


‘V - A A 

' A 




- 


X 

X X "X 

•A A A 

•A 

— 


- f v 

— 

- 

'•^ 

•v *s *v 

A *A 'A 

•A 


•S. 


V. 

V 

X 

v \ - V 

‘A V ‘A 

‘A 


**■ w. -V 


’X 

's 

V 

•\ V V 

A ’A A 

'A 

.t v 

■^v -X 

•X 

X 

V 

V 

■v -V V 

"V *A *\ 

'A 


•*X "X 

X 

V 

"V 

"V 

X X X 

X X X 

A 

^A 

\ V 

X 

X 

X 

X 

xxx 

X X X 

A 

% 

X X 

X 

X 

X 

X 

X X X 

X X \ 

\ 

\ 

X *\ 

\ 

X 

-r. 

\ 

■\ 

X -\ \ 

X \ V 

X 

\ 

\ X 

\ 

X 

X 

X 

X X X 

X X X 

X 

* 

\ \ 

\ 

\ 

\ 

\ 

■\ X 

X X X 

X 

\ 

\ \ 

\ 

\ 

\ 

\ 

\ \ \ 

X X X 

X 

\ 

\ \ 

\ 

\ 

\ 

' 

' \ A 

X x X 

X 

* 

1 

1 fc 

* 

\ 

\ 

V 

\ 

V v \ > 

X X X 

X 

1 

t \ 


\ 

\ 

•* 

A \ \ 

‘•t X X 

X 

N 

i 

\ \ 

\ 



> 

’* % t 

X X X 

X 

s 

t 

t ? 

t 

1 

V 

\ 

; 

t ^ ’t 

X X X 

X 

! 

; ; 

! 

1 

1 

> 

1 ^ ‘1 

'\ \ t 

X 

1 

; / 

1 

( 

t 

l 

\ i ’; 

\ \ X 

1 

l 

/ 7 

t 

V 

1 

< 


» * f 

7 7 7 

7 

» 

1 

./ / 

I 

t 

1 

V 

\ 

*i 

i* r r 

7 i 

f 

/■ 

] / 

l 

f 

s 

1 

% 

i 

p' t 4 

V V V 

t 

/ 

/ 7 

} 

t 

r 

) 

4 > ) 

f V V 

V 

/ 

/ / 

/ 

/ 

4 

) 

;- > / 

> > V 

V 

/ 

/ > 

4 

} 

? 

;■ 

;• ;■ ;■ 

V V 

V 

/ 

/ / 

/ 

} 

r 

> 

*t r > 

ir f V 

V 

_ a_ /_ r. 

_ ) _ 

_ } 

_ t . 

_ } _ 

_ 4 _ 7 _ < 

_ ■/ t . 4 _ 

_ > 


: O 
1 1 


nMU 

N/pHo 


-64 


-63 


co cn 













parison of OWI Model and HRD Winds During Hum 




UAt.M.tAKKU 




^ V ^\- 

i . if . T i f **v vwvvvwv ~v^ 

fr^N^^VvV^V^VH w '(V v^' A' V\^V \ v\ 
V^V^VV VWV'W \ * V\ Vv 

wwvwv f w\\*v\\\ v\ 

^WVVW\W\\\\ v\ 
WWVW\W\\V *w - V V 

X\X\X\‘VV\>X'V\V\> V\ -vv 

JotS\\\W\\W-vv w 

jNWv\x\\n\n\ V\\ \\\\ 


Sfe: va. 

>■,« \»V» V„V«V* V»V» 

BM. V.v. Y.V. wx VVXV VVVV>w\>kSN>v^^>^ 
BMV V.V. VX VWu\»VWW*fc 

■Mu VX V.V. VVV.VV. VV\VS.VVS^\1^^^ 
\A.\A.\A.V\AVv.\ .V. ^ 


ItI 


DOoX\^N\\\ \\ \ \\ \\ \\ 
\\\ V\ \\ 
J^NjvN \\N \\ \\ 

^ At®v*v »» v\\\ w w 


ikV VV^S,VVv.\.>uVi^*NW>i^^ 
Iv Vv.v wv\ .vwsw^ j awtfp yn 
Lv>.^\. VXS.N.\.Vv>>u> H » > » <m ■ ■ « « 

lw&t*'t»'nV- s t.N*'t> s v.> ,»<» K v 


spi 

W/Avau 


*m 4 j^jsa/aaa WKKKKK 1 ! P 1* 1 

>r^y‘^' , ^^A^AiWAl\\\\\\ \ ti lt I > 1 
tt'A'AfA'Af/} u H m i 

^^aaaaaaaa aaaa tUUV,X\\\\\\ v mi 
/wyvy/y//V // At At AAA A A f f rrf }} [( ? i 1 > i i 

^^aaaaaaaaaaaaaaaaaaaaaa i/limtiln 
^^^aaaaaa/a aaaa aaaa ajiau ;m n m m 

*w^ »>VV AAA/A At A A A A A A A A A A A A A } At A })) \ I } 
.■^■sassaaaaaaaaaa a/a tftnn n >> >>> m m 
.•A v>v/.v;/;/y aa a aaaa aaaaaaa n > n n n >>; 


! \j*vV\\\Vm 
\®vfK\\\\\\ 
S»i£0\\\\v\ 

F,iJ 5 J JJ JJU\ 

L, j lUini 
IlIJ Jl 3 ' tl 
1 n 3 m m n n 
Hj ummi 
ill Minin 
ill nm \ itt 
ill 1 h tun it 

:: mum 

(m h tUiljj 


^s;ss ss sss ts snn n>> > >? n >nn n njHrtn 
jjsss ss sss ts??n?n}> >}?}>)}) >},* tiittkri > t 

SS/ ' S S / f / r / r } f } } } ) > } } } } V}if kfj| H ftif 

ssssssssssssn n }>> >> a f \trff} 1> t V 
L^^ww;;/ s/s ss // / n n >> n n 0>r>t ► m > > b 


'////s'//sV/>V//n /fj 0 j£j/ It jx> r> tt 


<?n in 




















r 'k :Vr 






lO (T) 







r\l 


27 comparison of OVVI Model and HRD Winds During Hurricane Floyd 1999 1 40 


J'dLia 


•' 

> 

J 

/ 

* * * j 

ft ft s: : / A A 

A /■■■ A, a: A 

A 

A A s *> '. A.\ >' A A. A A — 

^ „ _ 

_ 

_ 



- 

— 


■ 

1 ^ 


' 

' 

/ 

ft ft / 

A / / * A A- A 

A A* s A M 

A 

Ak A. ft*-' A1 — 

- — ^ - 


• 

- 




_ 

a 

-i - 



' 

' 

/ 1 / /' 

A / / 

/ / '• A A A 

-* 

** 

A A A A ✓; <r A. A . **' A: 



- 

- 

- 


- 

- 

a 



' 

' 

/! 

/: / /■;: /■ 

/ / / A , / 

si /i & M M 


Aft Aft Ai A* Az A ^.; < i ** ■ A* >j . 

; — - — 

- 


- 

- 

- 

- 


a 

H 1 

/ 

' 

/ 

/ 

/ / / r 

* / / / / / 

A A A A A 

A 

Ai At Aft rfii ^ #<r. . At ft**-.- ^ 


- 

- 

- 

- 

- 

- 

- 



- / 

/ 

/ = 

/ ; 

' /i # /J 

A A / A A / 

A A A a a 

M 

iU Jiii >:•: > j riri -5> 



- 

- 

- 





5 

/ 

/ 

/; 


< v- £. a 

■ a a a ■ a 

a p ■ a ■ 



: a? -4. 

— 

— 

— 

— 

— 

• 

- 


■■3 

- 4 

/ 

/- 

/i 

/ ■; n a P 

• Mill 

ima 




— 

— 

— 

— 

^ . 

- 




1 / 

✓ 

> : 

/ 

/■ a ■ a 

Him 

aim 


M m £i M aA. A.- jt : .-r.' •— , •— •— 



- 

- 

- 

- 

- 

- 

- 

a 


L/ 

4 


a 

■ ■ * a 

Mini 

urn 


Jl ^ A)- A !*% r~.'- A r~ jml 



— 

- 

- 

- 

- 

- 

- 

a 


< 

4 

« 

a 

mama 

mm 

ppppp 


ji? 5«li •f‘5 ««*+' .*Ts, «!.: ^fc! ; 




-s. 





■ 

. '”3 

» * 

4 

■ 

a 

aaam 

aaaaaa 

ppppp 



4—: .. 


-S. 






am 

.5 


u> 


ft 




aaaaamaa 

aaaaamaa 






► '• N V \ 

niMliliiiii 

■ iiiMiiaa 


P p p a a. a. : 

■ ■■■■■»< 

■ ■■ 

P P P -«• - -- V 

■ P 

i|ll P '■ a *+-&' • Hr 'a 

i 'I *■*■••* ■ 

s *i *7 3 T 




/ j? -g t f: 

sr’ < 5 /’ V j i 

' i ii iimai 

1 A a a a a i i i ■ ■ 
aapaaaaiai 

tiviafiiiiimii 

1 :-r' ■> < > > i * 'iriiiiiii 

- * * *’ ,■? > if > > 


{ \ 

\ 

\ \ 

am 


aaaa 

aaam 

aaaa 


P P-P-PlP-P-P-Pa-PP-P-P-P * * 


HRD 

Mode 


/b 





























Comparison of OVVI Model and HRD Winds During Hurricane Floyc 

*i }.4 ’/ ^'/^^■/vrT'rrrry^rrrrr^ 


1999 


DSltiE 


" -« •< ^ y 

•» •» •■'■ ■' ' •' < * *f si&S «. **■ 

1 •* 1 it «6*J <£Sjjl£l |^£S +$&£ #££; 4H v ■*■*«+ + * + + + <+•*'*'«'■ 

•* ■* ■* „ im Iffi SSI tSI 2H£ ' !> " 

■t I 4 n ,cMb «W! MW ' 


. 

l 4 -M vm £M M K *p ^js , 




ISl^fiPP^lfliiU^ 


jfi^r’a >-.> , > 

J ; ji'i'! y y y / / jfl /aS S$ gfejlS ifiS? 

*•*«' y y y y y y y liejMfM 

yjy'yj, / y - y / , / ^W 4 aOT#il 

i.i j > > y y v j . . araa -*$£»$ Sis 


*~*-*~*-*— *-*-*-*-■•-*- -*- **-*•-*»- •- »* *« 
f'^'-f'p'rT-s*'* 3 -*-*'’**-*-*-*-*-*-*-**- 

r r T' t^f" F'P' fc *■*■*• *^ •>■ *■ 

-r^r- *“»*»• '*•*■*■ «*. 

r^r*- r-s>->- *»-*»■ ^ ^ ^ v 

^pjp- *~r*;z r r^” r=»t=“ j=- ?=* ft* *** ^ *« *« v 

*v Jk W 


^ ^ ^ A. ^ 

' ' - - - » ^ 

{jjiki*, -F ix ix -ix. ^ *X. 
..., - r .%V*u F, a. -n ^ 

-Zi»y^^3Nj£llfe^xAx AAA A A 

-■Ay.A-^U.^2 A A A A A A, A 
k«A*'\-~uA. A A AAA A A, 

k^\*AA. A A A A A A A A 

A*AAA A A A A A A A A A 
xAAAA A A AAA A A 


•VS -VV V - 




-1KJ 

Mcael 


-82 


-/o 


(jO Ul 



Comparison. g.f oyyi Mods! and HRD Wind^ n 



1 3 
11 


-15 

HRD 

Model 


-72 

! ^ 


-82 


-80 




-80 

















■■■■■■■■■a 


pmparison of QWI Model and HRD Winds During Hurricane Floyd 1 999 


4 4 ■* 4 4 :. 

4 : 4 4 4 4 :.. 


4 4 4 4 4 
4 4 4 4 4 

: 4* & /■ 4~ 4 l 

4 4 4 4 4 

4 4 4 4 4 
4 4 1 A '. A 4 
A' A 4 A A 


A A A A A 
A A A A A 
A A A \ A 
A A A A ( \ 
A: A A A A 
A A .% A 

A A ; \ 


/ y 4 s s s . . 

y / / 4 / / . ' s 

4 4 4 4 4 4 4 . 

y J 4 4 4 / 4 4 

J J J 4 4 4 X 

/ / J 4 , / . - 

/ / / / , 

l J j / / y / / . 







* vf '.“XU^ -V - 




j-V*\ ■ 

*A. * V* Vf <f - 


raaaa; 

r ■■■**.; 


*ia« 

.MBs'- 

T) ■ 


' {\ r V / 


V V Vll M 
< 1 V ( < **l 


N ^ ^aaa B wr >^ 1 1 t J J |ii :} :« j 

l*>^S5asaS5 


■ 'J •■. r w m T-T-T-T.*** , * «^a **! ^li ,t.i 

} fr4^vU.u.u .i XfeX. 1 ■* '} 
¥ M'44 >K'*VU •< - .} >.iu .* ‘. l 
;-Hv4'>AAA. -14 .1 ..' .J ..' ! 


. am<w r- 


«aaaaaa 

iiiassi 


- / / / / a 

'>*////// 


l«W 


AAA.i.i.i.t j .hi;::; 

V H < 'I I I 1 1 1 A i i 

® H *! * ** ** 7 7 7 

M ^ V 7 7 7 7 “i *t a 

w / / ** ,j * ** '* ^ ~ * 7 

■ i / / ^ V 7 7 7 7 V 7 

8 / J* ? V V 7 7 V V 7 7 

vti? y y y v v v v 7 

■ ■ S V 7 V V V V V V V 

J J r V V V V 7 7 7 V V 

g| ^ i' V V V V 7 V V V V 

7 / / / / / V V V V V V 


•sss////}}}} t / ^ V V V 


* * > d 


^ S ^ S S 4 / / / 
~ ^ S’ S’ 4 / / / / 


4 4 /> 4 t 4 


■■■■■■ 

vnaaiia- 


u ^ »■•■!?? : : ; ssshs 

4 s 4 4 / s / / / / ; > TT > f K ^ f V > V 2 S S 5 H 
m 4 4 44 ///t>i r ¥ k r ¥ t r r r V V Aj ^ ^ tljtl M 
* 4444 /}>t ¥ > r r r t r r V r ^ r M ti t tl Mi 

ttt*>¥>>¥¥VV \ ' llllllll 

A 4 4 4 > > k ¥ ¥ ¥ ¥ ', ¥ 

4 4 4 4 4 4 4 4 ¥ ¥ ¥ ¥ ¥ ¥ > 

a*44444¥¥¥¥>¥¥¥ \ ^ Mi ft ft H £ 
A/44444¥¥¥¥¥¥¥ M M ft ■ ■ M M ' 


-80 










Appendix D. Abstract of paper to 2001 AMS Annual meeting: 


Real Time Forecasting system of Winds, Waves and Surge in Tropical Cyclones: A Case Study 
with Hurricane Georges. 



Real-Time Forecasting System of Winds, 

Waves and Surge in Tropical Cyclones: 

A Case Study with Hurricane Georges 

Hans C. Graber, Mark A. Doneian and Michael G. Brown 

Rosenstiel School of Marine and Atmospheric Science 
University of Miami 

Peter G. Black , Mark D. Powell and Sam H. Houston 
NOAA/AOML Hurricane Research Division 
Vincent J. Cardone, and Andrew T. Cox 
Oceanweather, Inc. 

Robert E. Jensen 

• US Army Corps of Engineers 

Donald N. Slinn 

Florida Atlantic University 
John L. Guiney, 

National Hurricane Center 
Christopher Paluszek 
Sun Microsystems, Inc. 

The long-term goal of this partnership is to establish an operational forecasting system of winds, 
waves and surge impacting the coastline during the approach and landfall of tropical cyclones. 
The results of this forecasting system would provide real-time information to the National 
Hurricane Center during the tropical cyclone season in the Atlantic for establishing improved 
advisories for the general public and federal agencies including military and civil emergency 
response teams. As more people and societal infrastructure concentrate along coastal areas, the 
United States is becoming more vulnerable to the impact of tropical cyclones. Furthermore, it is 
not surprising that hurricanes are the costliest natural disasters because of the changes in the 
population and the national wealth density or revenue. A better understanding of both hurricane 
frequencies and intensities as they vary from year to year and their relation to changes in 
damages is of great interest to scientists, public and private decision makers and the general 
public. 

The estimation of tropical cyclone-generated waves and surge in the coastal waters and 
nearshore zone is of critical importance to the timely evacuation of coastal residents, and the 
assessment of damage to coastal property in the event that a storm makes landfall. The model 
predictions of waves and storm surge in coastal waters are functionally related and both depend 
on the reliability of the atmospheric forcing. Hurricane Georges is an excellent example of an 
intense tropical cyclone with numerous landfalls and unexpected changes in intensity and 
movement. Results of real-time forecasting simulations for Hurricane Georges will be presented 
as the storm segment moved over the Florida Straits and then entered the Gulf of Mexico. As it 
approached the coast of the central Gulf of Mexico it gradually slowed. Georges made its final 
landfall near Biloxi, Mississippi early on 28 September with 105 mph winds. Results will include 
both the predicted wave heights and surge-induced water levels along the Mississippi/Alabama 
coastline prior and during landfall of Georges. We will also show the variability of these 
parameters for different forecasted tracks and how such information would impact the advisories. 
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